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ABSTRACT 

We consider several protostellar systems where either a precessing jet or at least two misaligned jets 
have been observed. We assume that the precession of jets is caused by tidal interactions in noncoplanar 
binary systems. For Cep E, V1331 Cyg and RNO 15-FIR the inferred orbital separations and disk 
radii are in the range 4-160 AU and 1-80 AU, respectively, consistent with those expected for pre-main 
sequence stars. Furthermore, we assume or use the fact that the source of misaligned outflows is a 
binary, and evaluate the lengthscale over which the jets should precess as a result of tidal interactions. 
For T Tau, HH1 VLA 1 /2 and HH 24 SVS63, it may be possible to detect a bending of the jets rather than 
'wiggling'. In HH 111 IRS and L1551 IRS5, 'wiggling' may be detected on the current observed scale. 
Our results are consistent with the existence of noncoplanar binary systems in which tidal interactions 
induce jets to precess. 

Subject headings: accretion, accretion disks - binaries: general - stars: pre-main sequence - 
ISM: jets and outflows 
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been observed and calculate the parameters of the binaries 
in which tidal interactions would produce the observation- 
ally inferred precession frequencies. We next study cases 
where misaligned jets have been observed and, assuming 
or using the fact that the source of these outflows is a bi- 
nary, we calculate the precession frequency that would be 
induced by tidal interactions in the binary and the length- 
scale over which the jets should 'wiggle' or bend as a result 
of this precessional motion. In § ^ we give a summary and 
discussion of our results. 

2. THEORY OF PRECESSING WARPED DISKS 

We consider a binary system in which the primary and 
the secondary have a mass M p and M s , respectively. We 
suppose that the primary is surrounded by a disk of radius 
R with negligible mass so that precession of the orbital 
plane can be neglected. The binary orbit is assumed cir- 
cular with radius D and is in a plane with inclination angle 
5 to the plane of this disk. In general, S will be evolving 
with time. However, this evolution, which is not necessar- 
ily toward coplanarity, occurs on a long timescale (equal 
or larger than the disk viscous timescale) if the warp is not 
too severe (Papaloizou & Terquem 1995), so that we will 
consider here 6 as a constant. This means that although 
the outer parts of the disk may be driven out of the ini- 
tial disk plane on a relatively short timescale, the inner 
parts will retain their orientation with respect to the or- 
bital plane for a timescale comparable to the disk viscous 
timescale. Since jets are expected to originate from the 
disk inner parts, their orientation relative to the orbital 
plane will be determined by that of the disk inner regions. 

The secular perturbation caused by the companion leads 
to the precession of the disk about the orbital axis, as in 
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1. INTRODUCTION 

Most T Tauri stars are observed to be in multiple sys- 
tems (Mathieu 1994). There is also indirect evidence that 
the sizes of circumstellar disks contained within binary 
systems are correlated with the binary separation (Oster- 
loh & Beckwith 1995, Jensen, Mathieu & Fuller 1996). 
This suggests that binary companions are responsible for 
limiting the sizes of the discs through tidal truncation (Pa- 
paloizou & Pringle 1977, Paczyhski 1977). There are in- 
dications that, in binaries, the plane of at least one of the 
circumstellar discs and that of the orbit may not necessar- 
ily be aligned. The most striking evidence for such non- 
coplanarity is given by HST and adaptive optics images 
of HK Tau (Stapelfeldt et al. 1998, Koresko 1998). Also, 
observations of molecular outflows in star forming regions 
commonly show several jets of different orientations ema- 
nating from an unresolved region the extent of which can 
be as small as - 100 AU (Davis, Mundt & Eisloffel 1994). 
These jets are usually believed to originate from a binary 
in which circumstellar disks are misaligned, and in some of 
these systems a binary has indeed been resolved. In such 
cases, we expect tidal interaction to induce the precession 
of the disk (possibly both disks) which is not in the orbital 
plane, and thus of any jet it drives. Furthermore, a num- 
ber of jets seem to be precessing (see below). The above 
discussion leads us to interpret this precessional motion as 
being driven by tidal interactions between the disk from 
which the jet originates and a companion on an inclined 
orbit. 

In this Letter we consider several observed systems in 
the light of this model. In § we review the theory of 
precessing warped disks. In § 13 we apply it to observed 
systems. We first consider cases where a precessing jet has 



2 



a gyroscope. The disk is expected to precess as a rigid 
body if it can communicate with itself through some phys- 
ical process on a timescale less than the precession period. 
In non-self gravitating protostellar disks, communication 
is governed by bending waves (Papaloizou & Lin 1995, 
Larwood et al. 1996, Terquem 1998). The condition for 
rigid body precession is then satisfied if H/r > \u) p \ /fi , 
where u p is the (uniform) precession frequency in the disk 
and Slo is the angular velocity at the disk outer edge (Pa- 
paloizou & Terquem 1995). An expression for oj p has been 
derived by Papaloizou & Terquem (1995). Here we just 
give an approximate expression, that we derive by assum- 
ing that the disk surface density is uniform and that the 
rotation is Keplerian: 

15 M s (R\ 3 . IGM P 

where G is the gravitational constant. The assumptions 
under which this expression is valid (see Terquem 1998 for 
a summary) are usually satisfied in the case of the rela- 
tively wide binaries we will be discussing here and for the 
values of D/R we will be considering. Although we have 
assumed a circular orbit, an eccentric binary orbit can be 
considered by replacing D by the semi-major axis and mul- 
tiplying the precession frequency by (1 — e 2 )" 3 / 2 , where e is 
the eccentricity. On this basis we consider that the discus- 
sion presented below should remain valid for < e < 1/2. 

3. APPLICATION TO PARTICULAR SYSTEMS 

3.1. Expected parameters of binary systems with 
precessing jets 

Observations of molecular outflows in star forming re- 
gions show in some cases "wiggling" knots (or a helical 
pattern in projection onto the plane of the sky), which 
can be interpreted as being the result of the precession of 
the outflowing jet. In this section we will assume that such 
precession is caused by tidal interaction between the disk 
from which the outflow originates and a companion star 
in a noncoplanar orbit. In cases where the outflow has 
lasted many precession periods the implication is that the 
"wiggling" should be periodic. If we indeed assume that 
this is the case, the observations give the projected wave- 
length Xproj- When the angle i between the outflow and 
the line of sight can be estimated, the actual wavelength 
A = Xproj / sini can be derived. The precession period is 
then given by T = X/v, where v is the outflow velocity, and 
|w p | = 2it/T. Furthermore, if the outflow is precessing be- 
cause the disk plane and that of the orbit are misaligned, 
the angle 8 between these two planes is equal to the angle 
between the central flow axis and the line of maximum de- 
viation of the flow from this axis. This angle can also be 
observed. In all the cases studied in this section, 8 is small 
enough (~ 10 — 20°) so that we will consider cos 8 = 1. 

In this section we will adopt M p = 0.5 Mq. Since we 
do not know whether the jet which is observed to precess 
originates from the primary or the secondary, we will con- 
sider mass ratios M s /M p between 0.5 and 2. The lowest 
(largest) values would correspond to the case where the 
jet originates from the primary (secondary). These values 
are typical for pre-main sequence binaries. Finally, the 
disk from which the outflow originates is expected to have 
its size truncated by tidal interaction with the companion 



star in such a way that D/R lies between 2 and 4. Since 
Larwood et al. (1996) have shown that tidal truncation 
is only marginally affected by lack of coplanarity, in this 
section we will consider 2 < D/R < 4. Assuming a fixed 
ratio R/D, equation ([[]) can then be used to calculate the 
disk radius: 

"7*y. ( 2 , 

^32 M p \co p \ J \DJ v ' 

Since M p appears with the power 1/3, an uncertainty of 
a factor 2 over M p is equivalent to an uncertainty of only 
a factor 1.26 over R. The main uncertainty over R comes 
through the ratio R/D. We now consider some particular 
protostellar systems in the light of the above discussion. 
Cep E, at a distance of 730 pc, drives two outflows almost 
perpendicular to each other Eisloffel et al. (1996), which 
suggests that this source is a binary. In addition, they 
have interpreted the morphology of one of these jets as 
due to precession, and they have inferred T — 400 years. 
Figure [l].a shows a plot of D against R, as calculated from 
equation (0). We see that R lies in the range 1 — 10 AU 
while D lies in the range 4 — 20 AU. Binary separation 
would be 0."005 to 0."03, which is not currently resolv- 
able. 

V1331 Cyg is located at a distance of 550 pc. Visi- 
ble line emission shows a very faint and diffuse feature in 
the vicinity of this object, which appears to be a strongly 
'wiggling' jet (Mundt & Eisloffel 1998). The observations 
give X pro j — 0.5 pc (a full period is observed). Using 
i ~ 42° (Mundt & Eisloffel 1998), we derive A = 0.71 pc. 
Since v ~ 300 kms" 1 (Mundt & Eisloffel 1998), we get 
T = 2, 300 years. Figure [l|.b shows a plot of D against R, 
as calculated from equation (0). We see that R lies in the 
range 3 — 33 AU while D lies in the range 13 — 66 AU. Bi- 
nary separation would be 0."02 to 0."1. This upper value 
may be possible to resolve with the VLA or adaptive optics 
in the near-infrared. 

RNO 15— FIR, located at a distance of 350 pc, drives a 
molecular outflow which appears to be 'wiggling' (Davis et 
al. 1997). It is possible to interpret this morphology as due 
to precession within the uncertainty of the measurement 
(see Fig. 7 of Davis et al. 1997). From the observations, we 
derive X pro j — 0.065 pc. Assuming i — 45° (Cabrit 1989) 
and v = 10 kins" 1 (Davis et al. 1997), we get A = 0.092 pc 
and T = 9, 000 years. Figure [|.c shows a plot of D against 
i?, as calculated from equation (||). We see that R lies in 
the range 8 — 82 AU while D lies in the range 33 — 165 AU. 
Binary separation would be 0."09 to 0."47. Such a separa- 
tion may be possible to resolve with the VLA or adaptive 
optics in the near-infrared. 

3.2. Expected precession in binary systems with 
nonaligned jets 

In the systems presented in this section, misaligned "bi- 
nary" jets have been observed. Since it is very unlikely 
that one single source can drive two jets with very differ- 
ent orientations, it is assumed that each of the outflows 
originates from its own component of a binary system. In 
some cases a binary has actually been resolved, in other 
cases observations only allow us to put an upper limit on 
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the separation of the hypothetical binary. Since the out- 
flows arc not parallel, it is probable that the disks which 
surround these sources are themselves misaligned. There- 
fore, at least one of these disks is not in the orbital plane 
and should precess. We evaluate here the precession period 
T and give an estimate of the lcngthscalc A = v/T over 
which the outflows should 'wiggle' as a result of this pre- 
cessional motion. Since in general we do not know from 
which member of the binary each jet originates, we will 
here again consider mass ratios 0.5 < M s /M p < 2, unless 
otherwise specified. We will also take M p = 0.5 Mq un- 
less otherwise specified, and cos 6=1 (the results can be 
scaled for different values of 6 since lo p oc cos S) . 
T Tau is a binary located in Taurus, at a distance of 
140 pc. Observations show that two almost perpendicu- 
lar jets originate from this system (Bohm & Solf 1994). 
A disk of estimated radius R ~ 27-67 AU has been re- 
solved around the visible component, T Tau N (Akeson, 
Kocrncr & Jensen 1998). Here we assume that the disk 
around T Tau N is not in the orbital plane. We take 
D = 102 AU (Ghez et al. 1991) and for R the observa- 
tional values reported above. We fix M p = 0.7 Mq. Since 
we are interested in the precession of the jet emanating 
from the primary, we consider 0.5 < M s /M p < 1. We 
then get T ~ 5,000-4 x 10 4 years. Since v = 70 km s" 1 
for the jet emanating from T Tau N (Eisloffel & Mundt 
1998), A ~ 0.4-3 pc for this jet. These values of A are 
larger than the scale over which the jet has been observed 
so far (which about 0.1 pc), so that a bending rather than 
'wiggling' may be detectable in that case. 
HH1 VLA 1/2 is located at a distance of 480 pc. The 
two sources VLA 1 and 2, separated by 1,400 AU, drive 
the two misaligned jets HH 1-2 and HH 144, respectively 
(Rcipurth et al. 1993). We assume that this system is 
bound and noncoplanar, and that tidal truncation has op- 
erated such that 2 < D/R < 4 (note however that, in 
such a young and wide system, D/R may actually be 
significantly larger). We also assume that the angle be- 
tween the line of sight and the orbital plane is 45°, so 
that D ~ 1,980 AU. Then T ~ 4 x 10 5 -4 x 10 6 years. 
Since v ~ 200 km s" 1 for both flows (Eisloffel, Mundt & 
Bohm 1994), we get A ~ 77-870 pc. Since T is probably 
comparable to the age of the system, a bending rather than 
wiggling of the jet may be expected on a scale of a few pc. 
We note that 'wiggling' or bending of the jets has been 
suggested on the current observed scale, which is about 
0.5 pc for both jets (Reipurth et al. 1993). This clearly 
cannot be due to interaction between VLA 1 and 2, but it 
may be the sign that this system contains more sources. 
This is supported by the existence of at least two more 
outflows with different orientations (Eisloffel et al. 1994). 
HH 111 IRS: Perpendicular to the HH 111 jet, located 
at a distance of 480 pc in Orion, is another outflow called 
HH 121 (Gredel & Reipurth 1993; Davis et al. 1994). We 
assume that the source of these two outflows is a binary. 
From the unresolved central source in VLA observations 
(Rodriguez 1997) we infer an upper limit on the separa- 
tion D of about 0."1, or 48 AU. By adopting D = 48 AU 
and 2 < D/R < 4, we find that T ~ 1,000-8,000 years. 
Since v ~ 350 kms -1 for HH 111 (Reipurth, Raga & 
Heathcote 1992), A ~ 0.5- 6 pc for this jet. The low- 
est of these values is close to the extent over which the 
jet has been observed so far, which is 0.45 pc in pro- 



jection (Rcipurth 1989). We note that tidal effects in a 
putative binary system have been invoked by Gredel & 
Reipurth (1993) as a possible cause of the asymmetry of 
HH 121. 

HH 24 SVS63: The region of HH 24, located in Orion 
at a distance of 480 pc, contains several highly collimated 
outflows (Eisloffel & Mundt 1997) and a hierarchical sys- 
tem of four or even five young stars. The sources SSV 63E 
and SSV 63W are separated by about 4,500 AU in projec- 
tion onto the sky (Davis et al. 1997). SSV 63W is itself 
a binary. On the images taken by Davis et al. (1997) 
and given to us, we have measured that its projected sep- 
aration is 920 AU. We find that SSV 63E is probably a 
triple system: the projected separation is 350 AU between 
SSV 63E A and B and 975 AU between SSV 63E A and C. 
We will take these projected separations as indicative val- 
ues of the actual separations. At least two outflows with 
very different orientations originate from SSV 63E. These 
are HH 24 G (Mundt, Ray & Raga 1991) and HH 24 C/E 
(Solf 1987; Eisloffel & Mundt 1997) , which extend over 
0.2 and about 1 pc, respectively. SSV 63W is the source of 
another parsec-scale outflow, HH 24 J (Eisloffel & Mundt 
1997). Here again we fix 2 < D/R < 4. The velocities 
of the jets, in km s _1 , are about 140, 180 (if we assume 
the radial and tangential velocities to be similar) , 370 and 
50 for HH 24 J, HH 24 G, HH 24 C/E and for the red- 
shifted lobe HH 24 E, respectively (Jones et al. 1987). The 
different interactions within the two systems then lead to 
A between 5 and 556 pc, which is at least one order of 
magnitude larger than the extent over which the jets have 
been observed so far. This indicates that a bending rather 
than a 'wiggling' would be more likely to be observed. 
L1551 IRS5 is located in Taurus, at a distance of 
150 pc. It has been suggested that two jets could 
be emanating from this source (Moriarty-Schieven & 
Wannier 1991, Pound & Bally 1991). Furthermore, 
Rodriguez et al. (1998) have shown that this system is 
a binary with separation 45 AU and they have resolved 
two circumstellar discs for which they infer R ~ 10 AU. 
Their results are also consistent with the presence of two 
outflows which appear to be misaligned (see their Fig. 2). 
We take D between 45 and 63 AU, corresponding to an 
angle between the orbital plane and the line of sight in the 
range 0-45°, and we fix R = 10 AU. Then T ~ 4,000- 
5 x 10 4 years. Since for the jet which has been unam- 
biguously observed v ~ 200 km s _1 (Sarcander, Neckel & 
Elsasser 1985), we derive A ~ 1-10 pc. The projected ex- 
tent over which the jet has been observed so far is about 
1 pc (Moriarty-Schieven & Snell 1988). Therefore, even 
though a full period may not yet be seen, 'wiggling' could 
already appear in the observations. We note that the out- 
flow as observed appears to be very complex, and it may 
well be seen to precess. 

4. DISCUSSION AND SUMMARY 

In this Letter we have considered several protostellar 
systems where either a precessing jet or at least two mis- 
aligned jets have been observed. In the case where a jet 
is seen to precess (or rather interpreted as precessing), we 
have assumed it originates from a disk which is tidally 
perturbed by a companion on an inclined orbit, and we 
have evaluated the parameters of the binary system. For 
Cep E, V1331 Cyg and RNO 15-FIR, we found the sepa- 
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ration to range from a few AU up to 160 AU and the disk 
size to be between 1 and 80 AU. These numbers corre- 
spond to what is expected in pre-main sequence binaries 
(see Mathieu 1994 and references therein). We note that 
larger separations for this range of disk sizes would be 
associated with longer precession timescales, and thus it 
would be more difficult to detect the precessional motion 
over the observed lengthscale of the jet. A bending rather 
than 'wiggling' would be expected in that case. 

In the case where misaligned jets have been observed, 
we have assumed or used the fact that the source of these 
outflows is a binary, and we have calculated the precession 
frequency that would be induced by tidal interactions in 
the binary and the lengthscale over which the jets should 
'wiggle' as a result of this precessional motion. For T Tau, 
HH1 VLA 1/2 (which may actually be a hierarchical sys- 
tem) and HH 24 SVS63, it may be possible to detect a 



bending of the jets rather than 'wiggling' on the current 
observed scale of the jets. In HH 111 IRS and L1551 IRS5 
(assuming there are indeed two misaligned outflows in this 
system), 'wiggling' may be detected on the projected scale 
(0.5-1 pc) over which the jets have been observed. 

Our results are consistent with the existence of non- 
coplanar binary systems in which tidal interactions induce 
jets to precess. Some of the predictions of this Letter could 
be tested observationally in a near future. 
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on the Dynamics of Astrophysical Discs, when this work 
began. CT is supported by the Center for Star Formation 
Studies at NASA/ Ames Research Center and the Univer- 
sity of California at Berkeley and Santa Cruz. 
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Fig. 1. — The binary separation D vs. the disk radius R in AU, as calculated from equation (p2), for Cep E (a), V1331 Gyg (b) and 
RNO 15-FIR (c). The different curves correspond to M a /M p = 0.5 (solid line), 1 (dotted line) and 2 (dashed line). D/R is kept in the range 
2-4 and the symbols indicate D/R = 4 (crosses), 3 (squares) and 2 (triangles). 



